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Abstract 

In models with large extra dimensions, where the fundamental gravity scale can 
be in the electroweak range, gravitational elFects in particle physics may be noticeable 
even at relatively low energies. In this paper, we perform simple estimates of the decays 
of elementary particles with a black hole intermediate state. Since black holes are 
believed to violate global symmetries, particle decays can violate lepton and baryon 
numbers. Whereas previous literature has claimed incompatibility between these rates 
(e.g. p-decay) and existing experimental bounds, we find suppressed baryon and lepton- 
violating rates due to a new conjecture about the nature of the virtual black holes. 
We assume here that black holes lighter than the (effective) Planck mass must have 
zero electric and color charge and zero angular momentum - this statement is true 
in classical general relativity and we make the conjecture that it holds in quantum 
gravity as well. If true, the rates for proton-decay, neutron-antineutron oscillations, 
and lepton-violating rare decays are suppressed to below experimental bounds even 
for large extra dimensions with TeV-scale gravity. Neutron-antineutron oscillations 
and anomalous decays of muons, r-leptons, and K and i3-mesons open a promising 
possibility to observe TeV gravity effects with a minor increase of existing experimental 
accuracy. 
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1 Introduction 



In physics, we sometime rmi into the so-cahed hierarchy problems. This happens when we 
find two or several quantities which are different by many orders of magnitude whereas, a 
priori, we would have expected them to be more or less at the same scale. 

A very well known hierarchy problem in high energy physics is related to the huge gap 
between the Planck mass Mpi ~ 10^^ GeV, which sets the energy scale for gravitational 
interactions, and the electroweak scale Mew ~ 10'^ GeV, which instead is joined with 
the particle physics world. In recent years, models with extra dimensions have been pro- 
posed [T]-[4], where the true fundamental gravity scale can be as low as a few TeV and where 
the Planck mass is a mere effective long-distance 4-dimensional parameter. For a recent 
review see [5]. In these scenarios, gravity becomes phenomenologically interesting for high 
energy physics and we may have realistic possibilities to observe and study quantum gravity 
phenomena at future colliders. 

However, one of the constraints on theories with large extra dimensions is the rate of rare 
processes such as proton decay. Current experimental bounds depend on the decay mode 
in question [6], with several channels bounded by Tp > 10^^ years [71 18]. Contributions to 
p-decay p] include those from GUT-scale intermediate bosons that mediate baryon number 
violation; these contributions can be suppressed below experimental bounds by additional 
symmetries as discussed by [H [21 [lOl [TTl [12l [13] . In addition, several authors have worried 
that proton decay and other rare B- and L-violating decays via virtual black holes (BHs) can 
be exceedingly rapid. According to common belief, the decay/evaporation of BHs does not 
conserve any global C/(l)-quantum numbers and, in particular, baryonic, B, and leptonic, 
L, charges [H]. The usual description of this phenomenon is based on enumeration of all 
possible operators which do not conserve global charges, normalized to the Planck scale 
with dimensionless coefficients of order unity. 

The suggestion that proton decay can be mediated by a virtual BH was put forward 
about 30 years ago by Zeldovich [15] and subsequently considered by [MljTTl [13 [l9] . The 
probability of this process in the standard frameworks with Mpi ~ 10^^ GeV is very low, 
but with smaller Planck mass in TeV range, the baryon number violating processes through 
formation of an intermediate virtual BH become much more efficient. Adams et al. [20j 
argued that experimental limits on the proton lifetime constrain the quantum gravity scale 
to be larger than 10^^ GeV, implying that the size of the large extra dimensions is less 
than / < 10^/" X 10 cm, where n is the number of such dimensions. Arkani-Hamed and 
Schmaltz [21j propose reducing the proton decay rate by a scenario in which the Standard 
Model fields are confined to a thick wall in extra dimensions, with fermions "stuck" at differ- 
ent points in the wall; then couplings between them are suppressed due to the exponentially 
small overlaps of their wave functions. 

Here, instead, we suggest that p-decay via virtual BHs may be more suppressed than 
previously thought. We propose a "conjecture" that, just as in classical gravity, sub-Planck- 
mass BHs in quantum gravity can only exist with zero local charge (electric or color) and 
zero angular momentum. In classical general relativity, it can be shown that a BH with mass 
smaller than the (effective) Planck mass must have zero electric charge and zero angular 
momentum (see Section 4); otherwise no horizon can exist and the naked singularity is 
exposed. Moreover, it can be shown that it is impossible to make the electric charge of 
a classical BH larger than its mass, > M^, in Planck units, by an influx of charged 
particles on the horizon (see e.g. \22\ I23|). In this paper we make the (as yet unproven) 
conjecture that this result remains true in quantum gravity. Then the virtual BHs that 
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mediate interactions such as proton decay must have zero spin and be electricahy neutral. 
As a consequence, a four-body colhsion (3 quarks and a lepton) is required for formation of 
the BH intermediate state, which leads to a suppressed proton decay. 

We compute rare baryon and lepton violating decays as well as (n — n)-oscillations 
due to intermediate virtual BHs and find that these rates are in agreement with existing 
experimental bounds even for TeV scale fundamental gravity. In preparation is a second 
paper, in which we will investigate a mechanism for generating the baryon number of the 
universe using BH-mediated processes. 

The content of the paper is as follows. In Sec. [2] we briefly review TeV-scale gravity 
models and in Sec. [3] the standard approach to gravitationally induced proton decay. In 
Sec. m we discuss our BH conjecture. In Sec. [5] we consider rare lepton-violating decays 
with a BH intermediate state, giving an estimate of their rates. In Sec. [6]we discuss baryon- 
violating processes, proton decay and (n — n)-oscillations, with a BH intermediate state. In 
Sec. [7] some anomalous decays of particles with heavy quarks, t and b, are evaluated. The 
results are summarized in Sec. [HI The paper closes with an appendix where we review 
present lower bounds on the magnitude of the possible fundamental gravity scale M*. 

2 TeV-gravity models 

In 1998 Antoniadis, Arkani-Hamed, Dimopoulos and Dvali proposed a "geometric" solution 
to the hierarchy problem of high energy physics, where the observed weakness of gravity 
(at long distances) would be related to the presence of large compact extra dimensions 
[U |2] . Motivated by string theory, the observable universe would be a 4-dimensional brane 
embedded in a (4+n)-dimensional bulk, with the Standard Model particles confined to 
the brane, while gravity is allowed to propagate throughout the bulk. In such scenarios, 
the Planck mass Mpi becomes an effective long-distance 4-dimensional parameter and the 
relation with the fundamental gravity scale is given by 

Mj>i ~ M^+^R'', (1) 

where R is the size of the extra dimensions. If these extra dimensions are "large", i.e. 
R 3> Mp^ ~ 10~^^ cm, then the fundamental gravity scale can be as low as a few TeV and 
therefore of the same order of magnitude as Mew- If we assume M^^ ~ 1 TeV, we find: 

R io(30/")-i7 (2) 

In this approach, however, the hierarchy problem is not really solved but shifted instead 
from the hierarchy in energies to a hierarchy in the size of the extra dimensions which are 
much larger than 1/TeV ~ 10~^^ cm but much smaller than the 4-dimensional universe 
size. 

The case n = 1 is excluded because from Eq. ([2]) we would obtain R ~ 10^^ cm and 
therefore strong deviations from Newtonian gravity at solar system distances would result. 
For n > 2, R < 100 //m and nowadays we have no experimental evidence against a modi- 
fication of gravitational forces in such a regime [23]. Interesting variations of these models 
can lower the fundamental gravity scale with the use of non-compact extra dimensions . 

If gravitational interactions become strong at the TeV scale, quantum gravity phenom- 
ena are in the accessible range of future experiments in high energy physics. In particular, 
there is a fascinating possibility that hadron colliders (such as LHC) will be BH factories 
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(for a review, see e.g. Ref. [25) . criticisms can be found in Refs. |26t l2Uj). From the clas- 
sical point of view, we expect BH production in collision of two particles with center of 
mass energy ^/s if these particles approach each other so closely that they happen to be 
inside the event horizon of a BH with mass Mbh ~ ^/s- Semiclassical arguments, valid for 
Mbh ^ M^, predict the BH production cross-section 



ttRIjj{Mbh) 



(3) 



where Rbh{^bh) is the horizon radius of a BH of mass Mbh- In the case of an uncharged 
and non-rotating BH in (4+n)-dimensions, the horizon radius is obtained from the {A+n)- 
dimensional Schwarzschild metric [27J 
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where we have ignored possible effects of the gravitational field of the brane and have 
assumed Rbh much smaller than the size of the extra dimensions, so that the boundary 
conditions which come from compactification can be neglected. On the contrary, if Rs is 
larger than some extra dimensions, the Schwarzschild radius is given by a lower dimensional 
BH solution. We have introduced here two quantities, Rbh for the horizon radius of a 
generic BH and the same for a Schwarzschild BH. Since in what follows we will deal 
with Schwarzschild BHs only, these quantities are the same. 



3 Spacetime foam and proton decay 

In this section, we discuss the previous works [19^ [20] in which BHs from the spacetime 
foam give rise to proton decay. We raise some objections and, in the next section, turn to 
our conjecture which implies that the decays of this section do not take place. 

Since at the moment a reliable quantum theory of gravity is lacking, the standard 
approach is based on semiclassical calculations. For example, virtual BHs can be obtained 
considering the path integral for N non-interacting BHs which, in a 4-dimensional spacetime, 
is (the generalization to a (4+n)-dimensional spacetime can be found in Ref. [20] ) 

where y is a normalization volume and Lpi is the Planck length. From Eq. ([5]) we can define 
a probability distribution of having BHs with mass m and, computing the expectation 
value of their mass and of their number density, we find (m) ~ Mpi and a density of 
roughly one BH per Planck volume. This is the so called spacetime foam picture [301 El] 
where the spacetime is filled with tiny virtual BHs which arise as quantum fluctuations out 
of the vacuum. These BHs exist (for very short time, of order of one Planck time, as a 
consequence of the uncertainty principle) in the same sense as electron-positron pairs exist 
in the vacuum of quantum electrodynamics. 

In this context we consider the decay of the proton due to virtual BHs from the spacetime 
foam. We take the proton as an object of volume Vp equal to its Compton wavelength cubed, 

^In this paper we use the convention of Ref. [28|: M"'^^ = 1/G*, where G, is the (4+n)-dimensional 
gravitational constant. See e.g. [25] for other conventions. 
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that is Ip ~ Jn-p^, with three valence point-Uke quarks inside. Probably a better estimate is 
Vp ~ Ag^^, where Aqcd = 100 — 300 MeV. Thus in the estimates presented below Aqcd 
may be substituted instead of nip. The probability of finding two quarks in the same Planck 
volume is simply given by the ratio of the Planck volume to the box volume [nip/Mpi)'^ . If 
we include the fact that a virtual BH should be formed at the same time when these two 
quarks are packed so closely, we obtain an additional nip/Mpi suppression factor and, by 
dimensional analysis, we arrive at the decay rate 

Since BH evaporation conserves energy, charge and angular momentum, this expression is 
an estimate of the rate of the following process: 

q + q^q + l, (7) 

where q (q) is a quark (anti-quark) and I a charged lepton. 

The same result can be obtained by different arguments. The rate of the process ([7]) 
can be estimated as: 

h/n = naBH = crBH\tp{0)\'^, (8) 

where n ~ is the number density of quarks inside the proton and asH is their inter- 
action cross-section through formation of a virtual BH. Since the interaction arises from 
a dimension six operator, the amplitude has a factor l/Mp^ and the cross section can be 
estimated as 

asH ~ ml/Mj^i. (9) 

Again, we obtain the result in Eq. ([6]). 

Inserting the standard Planck mass Mpi ~ 10^^ GeV into Eq. ([6]) we predict the proton 
lifetime of the order of 10^^ years and there is no problem with the current experimental 
bound. On the other hand, in models with large extra dimensions the fundamental gravity 
scale is and, replacing Mpi with M^, ~ 1 TeV, leads to quite short-lived proton with the 
lifetime Tp ~ 10^^^ s. Hence, in order to avoid contradiction with the present experimental 
constraints, Adams et al. [20J require > 10^^ GeV. 

However, some criticisms may be raised in this connection. In fact, even if in some 
cases vacuum fluctuations lead to observable phenomena, such as the Casimir force, they 
are doubtless not well understood when we consider quantum gravity effects related to 
formation of small virtual BHs. We may encounter problems, even with the standard Plank 
mass of 10^^ GeV, if we apply the same considerations to leptons, which are supposed 
to be elementary ones (not consisting of some "smaller" parts). In this case much faster 
gravitationally induced decays can be expected, because we do not need to wait for two 
constituent particles to be within the same Planck volume. Hence, in the arguments leading 
to Eq. ([6]), only the m/Mpi suppression factor remains. For example, the rate of /i ^ 67 
could be as large as F ~ rnj^/Mpi. This estimate is based on time uncertainty relation 
and gives rise to the decay branching ratio about 10~^, which contradicts the experimental 
bound by 6 orders of magnitude. 

However, these considerations are pretty inaccurate and cannot rigorously lead us to the 
conclusion that the spacetime foam picture with the standard Planck mass is incorrect. In 
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particular, the result presented above for the rate of the decay /x — > 67, which is inversely 
proportional to the first power of the Planck mass looks very suspicious, one should expect 
at least the second power, and possibly indicates a sickness of this kind of argument, because 
in any normal theory the probability of a process in the leading order must contain an even 
power of the coupling constant. This is true in the usual approach, according to which 
the decay — > 67 is described by the dimension five operator, F^^{ipa'^'^ip) /Mpi so that 
r ~ m'^/M^i. It does not contradict the experimental bound with the normal Planck mass, 
Mpi = 10^^ GeV, but still it puts a very strong lower limit on the effective Planck mass, 
> 1013 GeV. 

In fact practically all "natural" estimates with the low scale gravity lead to too high 
rates of processes with baryonic and/or leptonic number violation. One possible conclusion 
is that the gravity scale must be high, close to Mpi, or that there is some new exclusion 
principle which forbids an easy formation of virtual BH. In what follows we will explore 
the latter option, requiring zero charge and angular momentum of virtual BHs because 
non-zero charge and momentum are classically forbidden for very light ones with m < M*. 
Simultaneously here enters another important assumption that virtual BHs have masses of 
the order of the masses of the initial state (see the next section). 



4 A Classical Black Hole Conjecture 

In classical general relativity, a BH with mass smaller than the (effective) Planck mass 
cannot be formed if it has a non-zero electric charge or if it rotates. Here we will make the 
conjecture that this result remains true in quantum gravity. Hence the virtual BHs that 
mediate interactions such as proton decay must be electrically neutral and spinless. 

The classical condition that a BH less massive than the Planck mass must have vanishing 
angular momentum and electric charge can be seen as follows. We can examine the 4 
dimensional Kerr-Neumann or Reissner-Nordstrom solutions to see that a horizon does not 
exist if the BHs have electric charge or angular momentum. The position of the horizon is 
given by the expression [22l [23] 

where Q and J are respectively electric charge and angular momentum of BH. It is clear 
that in 4D space-time the horizon cannot be formed if 

In the absence of a horizon, there would be naked singularities which are not allowed. 
Indeed if condition (jlip is fulfilled, the Kerr-Newmann metric allows to reach the physical 
singular point r = from some large r in finite time without crossing any metric singularity. 
Probably BHs with such low masses cannot be created. See e.g. problem 17.17 in book [23] 
where it is shown that it is impossible to charge an already existing Reissner-Nordstrom BH 
up to Q > Mbh/Mpi by any physical process. So a naked singularity cannot be created. 

This result is true only in classical physics and quantum gravity effects may change 
the bound (llip . Moreover, the time-energy uncertainty relation could allow a virtual BH 
to have mass/energy larger than the mass of the initial particle. Still in the absence of 
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anything better we assume that only diagrams with uncharged and non-rotating virtual 
BHs can participate in light (m < Mpi) particle decays. 

As for the color charge, we expect a similar picture, though the situation is much less 
clear. For classically large BHs any color charge should be screened at the microscopic 
distances of the order of the inverse QCD scale, Rqcd = I/^qcDj so there are no colored 
hairs at large distances but they can be present at the distances smaller than Rqcd- The 
analysis of classical colored BH properties can be found in Ref. For microscopically 

small BH, considered here, the screening radius Rqcd is usually large in comparison with 
the radius of the BH. Thus for sub-Planck-mass BHs we can apply the same prescription 
as is done above for electrically charged BH and require colorless states, even if an analytic 
bound such as Eq. (jllh does not exist. On the other hand, weak SU{2) charges can be 
safely neglected, since the symmetry is broken via the Higgs mechanism and BHs do not 
manifest charges if they are related to massive gauge bosons [33]. To be more precise, weak 
hairs may exist but only for a short time, tw ~ l/A^VKj where M\y is the mass of weak 
gauge boson, W or Z. Since this time is much shorter than the life-time of the processes 
considered below we neglect weak hairs of BHs and allow for arbitrary quantum numbers 
with respect to weak SU{2). 

A similar situation should hold in higher dimensions but instead of Mpi an effective 
gravitational scale should be substituted. If the Compton wavelength of an elementary 
particle, Ac = 1/m, is much smaller than the size of the extra dimensions, it is natural to 
believe that gravity "inside" an elementary particle becomes multidimensional. 

In addition to this conjecture of neutral and non-rotating BHs we impose some, maybe 
even more questionable, conditions in calculations/estimates of the amplitude of reactions 
with broken global quantum numbers due to virtual BH. In essence we suggest a set of 
rules which do not respect some of the usual conditions existing in quantum field theory, 
in particular crossing relations between amplitudes. For example, we allow a virtual BH to 
decay into, say, a proton and a electron, but we do not allow a proton to form a BH plus 
a positron, with the same amplitude. The picture that we have in mind is a kind of time 
ordering: a BH could be formed in a collision of a neutral system of particles in the s-channel 
whereas a BH cannot be in the t-channel of a reaction. We assume that BHs can be formed 
out of positive energies of real particles only and not from virtual energies of particles in 
closed loops. For example, BH cannot be formed by vacuum fluctuations, despite the fact 
that, according to the standard picture, vacuum fluctuations might create a pair or more of 
virtual particles both with positive and negative energies. The mass of the BH should be of 
the order of the energy of incoming (or outgoing) particles. In an attempt to describe this 
in terms of the usual language we come to a version of the old non-covariant perturbation 
theory with all virtual particles having positive energies. It corresponds to the choice of 
only one mass-shell pole in the Feynman Green's functions. This rule allows only for BHs 
with masses which are of the order of the energies of the initial (or final) particles, as we 
postulated above. It may look very strange, to say the least, but virtual BHs are not well 
defined objects and we do not know what happens with space-time at the relevant scales. 
Taken literally these rules would lead to violation of some sacred principles of the standard 
theory (locality, Lorentz invariance, and more). Let us remind the reader, however, that the 
existing attempts in the literature to invoke virtual BHs are based on standard quantum 
field theory in a situation where it is almost surely inapplicable. 

So it is not excluded that many properties of the standard field theory are broken, 
including even Lorentz invariance and locality. We cannot of course present any serious 
arguments in favor of our construction but it predicts quite impressive phenomena with 
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clear signatures based on a very simple set of rules and if these effects are discovered, the 
approach may be taken more seriously. Our goal here is to formulate a reasonable(?) set 
of rules which may possibly describe processes with virtual BHs and are, at least, not self- 
contradictory. Based on these rules we will study the phenomenological consequences which 
are quite rich and may be accessible to experiments after a minor increase of accuracy. 

5 Lepton number violating decays mediated by black holes 

Gravitational interactions are known to become stronger at short distances or at high ener- 
gies and in the standard theory they are expected to become strong at ~ Mpi ~ 10^^ GeV 
or at the distances about 10~^^ cm. In TeV-gravity models, we can expect non-negligible 
gravitational effects much earlier at the energy scales of contemporary particle physics, at 
relatively low energies, about electroweak scale. Let us remind the reader that rare de- 
cays have been often used to probe small distances and sometimes have given information 
about heavy particles prior to their discovery. They can play the same important role to 
investigate models with the fundamental gravity scale in the TeV range. 

5.1 Muon decay /i — > 3e 

In this connection we will consider /i~ decay into e~e~e^, which violates muonic/electronic 
number conservation due to formation of a virtual BH in the intermediate state. The 
diagram describing this process is presented in Figure [Hi. 

First, the muon emits a virtual photon; then the photon produces an electron-positron 
pair. Next the muon and the positron form a Schwarzschild BH. Since the BH does not 
respect muon-number conservation, it can decay into an e"^e~-pair. This is not the Hawking 
radiation |34j because the latter is a semiclassical process, which can be realized for "large" 
BHs, but is essentially a quantum gravity phenomenon. We cannot reliably calculate its 
amplitude, since it is surely non-perturbative, but assume that it can be estimated on 
dimensional grounds assuming that numerical coefficients are of order unity. We also assume 
that the BH decays predominantly into Standard Model particles on our brane [35] (however, 
since bulk emission of gravitons becomes more relevant as the number of extra dimensions 
increases ^36j, even processes with missing energy may play an interesting role). BH decay 
conserves energy, angular momentum and gauge charges, but violates global symmetries. 
As we have already mentioned the process presented in Figure [Th violates the conservation 
of lepton family number. 

Let us now make a rough estimate of the rate of this decay. The emission of the virtual 
photon and its subsequent transformation into an e+e~-pair leads to the suppression of the 
probability by the factor of the order of (a/27r)^. 

By dimensional arguments, the amplitude of the decay fj,^e^ into e^e~ through a virtual 
BH is proportional to ^I/^bh' where Mbh is the mass of BH and §2 is the dimensionless 
coupling constant of BH to two fermions. This coupling constant must be proportional 
to the strength of gravitational coupling. Here we make the assumption that §2 ~ RsE-, 
where E is the total energy of all colliding particles which make the virtual BH in their 
center of mass system, i.e. E = Mbh- Thus g2/MBH = Rs- As is written above, Rs is the 
Schwarzschild gravitational radius. For Rs we use Eq. ([4]) with the last factor in square 
brackets taken to be ^/^^ in the case of multidimensional gravity, while Rs = Mbh /Mp^ for 
the standard gravity. 
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We have identified the energy of the colhding or outgoing particles with the BH mass, 
while Rs describes the strength (or better to say, weakness) of gravitational interaction. 
We emphasize that this is an assumption about the nature of the (quantum) gravitational 
interaction with the BH. As for the mass of the BH we assume that it is of the order of 
the muon mass, M^h ~ t^^i^ or niore precisely, M^h = + E^- in their rest frame. In 
principle one might worry that a virtual BH might have a much larger mass, e.g. larger 
than M*, and in this case it could be both charged and rotating. For very heavy BHs in the 
intermediate state the amplitude should be suppressed by an inverse power of BH mass and 
hence it could be weaker than the amplitude with light neutral BH. On the other hand, if 
a heavy charged and/or rotating BH could be formed in low-body (e.g. two-body) collision 
of the constituents of the initial particle, while neutral BH demands more particles for its 
creation, the processes with heavy BHs would be less suppressed and could be dangerously 
efficient, despite an absence of the factor (a/27r)^. A much larger M^,, beyond TeV range 
would then be necessary in this case to avoid conflict with experiment. 

Quantum field theory allows for any masses of virtual particles. However, it may be 
not true for virtual BHs which are surely nonperturbative and quite complicated quantum 
fluctuations of space-time. We suggest that such heavy intermediate black holes are not 
allowed. Our hypothesis is that virtual BH's are in some sense real, i.e. their mass is of the 
order of the real characteristic energy of colliding particles which form BH in the process 
under scrutiny. This assumption contradicts the statement of quantization of BH masses 
with steps of the order of Af* but in the absence of a serious theory it may be as good as 
any other. At least, it is very simple and in some sense natural: for formation of BH we 
need real available energy which is of the order of the decaying particle mass. 

The amplitude of /i — > 3e decay corresponding to the diagram in Fig. la is equal to: 

a In (M^/ml) r 2 _ _ 2 _ _ i 

"^^^ ^ ^ vrM2 [(723)^ (^2^3)(^i^m) - (713)^ (V'1^3)(^2V'm)J , (12) 

where 'jij = M^^^^/M^ come from the factor gl/Ml^ = (see Eq. (g])) and a = 1/137; 

the log-factor comes from logarithmically divergent triangle part of the diagram with the 

ultraviolet cut-off taken at the gravity scale A.uv ~ V'j £^re the Dirac spinor wave 

functions of the corresponding particles. The amplitude is antisymmetric with respect to 

interchange of two final state electrons 1 and 2. The mass of the virtual BH is taken to be 

equal to the energy of the e~^e~ pair emitted by this BH. The upper indices -Af^^ indicate 

(ii) 

which particles are emitted by BH. In what follows we substitute for simplicity for M^jj 
its average value 

((^BH )') = (K - 2m^^3)) = Kml, (13) 

where k = 1/2 — 1/3. Taking the necessary traces and integrating over three body phase 
space we find the decay width: 

r(.-3e). = ^(ln-|j (^^j ..... (14) 

Two factors 1/2 come from averaging over the spin of the initial muon and from 1/2! due 
to two identical particles in the final state. 
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In the case of the ordinary (3+l)-dimensional gravity with = Mpi, the decay rate 
is roughly equal to 

n,^3eU,-0mi{^y^^ (15) 
and is negligibly small. In the higher dimensional case with n = 2: 

r(;. ^ 3e), - 6 • lO'^^ GeV (l + 0.11 In ' (^^^ ^ (3.)^/^ . (16) 

Since the total decay rate of the muon is 

Ftot « 3 • 10-19 GeV (17) 

then, assuming M* ~ 1 TeV, we obtain the branching ratio about 2 • 10^ for n = 2 which 
is close to the present experimental constraint [8J: 

<1.0-10-i^ (18) 

Exp 

For larger n the decay rate would be in stronger disagreement with the experimental bound. 
However even a minor increase in the value of M^, would avoid the contradiction. It is 
intriguing that these estimates with ~ TeV are quite close to the existing bounds on 
muon number violating decay /x 3e. 



BR{fi- e'e+e") 



5.2 Lepton number violation in e~^e collisions 

The conservation of muon number can be also violated in the reaction 

e+ + e"^/i + e (19) 

or similar reactions with any other leptons in the final state. A virtual BH can be formed 
here just from the initial e+e" pair and the cross-section of this reaction would be about 

^ ^ VlOOGeV J J ^ ^ 

This may be observable in high energy e'^C collisions. We again note that the predicted 
cross-section (I20p is surprisingly close to the current limits on lepton flavor violation search 
at the Z° resonance [37j . One should keep in mind, however, that the total angular mo- 
mentum of the initial e~^e~ pair must vanish. This means that the annihilation should 
proceed from the s-state because the total spin of the initial particles is zero. This demands 
the same sign of helicity of and e~ in their center of mass frame. On the other hand, 
the annihilation through Z-boson proceeds with total angular momentum equal one and 
demands opposite helicities of e"*" and e~. 



5.3 Muon decay ^ 

One can obtain an estimate of the decay rate /i — > 67 along the same lines. In the simplest 
diagram (see Figure 16) the probability is suppressed by an additional power of a, but the 
ratio of two-body to three-body phase space compensates this extra suppression, so that the 
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probability of /i ^ 67 decay through the considered mechanism would be approximately 
the same as — > 3e. 

We should keep in mind however, that there is an additional ambiguity related to the 
second loop in diagram of Fig. which contains three virtual electrons. If we first inte- 
grate over the {fie^)-loop, the second loop would also be logarithmic, because one electron 
propagator, common for both loops, would be reduced to a point. The log-factor related 
to the second loop should be taken with the loop counting factor, l/(27r)^ and the result 
would be somewhat smaller than the naive estimate without logarithmic and loop counting 
factors. There is an important difference between the first and the second loops. The par- 
ticles in the first loop interact with the BH at one point and the BH propagator does not 
enter into consideration. In the second loop the interaction with the BH is not pointlike 
and the statement about its ultraviolet behavior was made assuming the standard form of 
the BH propagator, ~ {Pbh ~ ■^bh)~^ ■ unknown whether or not this is true. Possibly 
covariant perturbation theory is not applicable to the case of virtual BH and naive estimates 
of divergences are not valid (see SecH]). This problem is also encountered in Sec. 16.31 where 
it is essential for the calculation of the time of (n — n)-oscillations. 

5.4 T decays 

The widths of the analogous processes with nonconservation of the r-flavor number in r- 
lepton decays would be enhanced by the ratio (rrir / m^)^~^^^ ^"^^^^ which is 6 • 10^ for n = 2. 
Since the r life-time is lO"^ times shorter than that of the muon, the expected branching 
ratios of the decays r — > 3Z or r ^ ^7 would be an order of magnitude larger than those 
for similar muon decays, i.e. maximally around 10^^^. It surely does not contradict the 
existing bounds, which are of order of 10~^ — 10~^ [8]. It is unclear if such decays may be 
reachable in not too distant future. 

If the virtual BH in the diagram in Fig. (TJi emits a gg-pair, instead of a lepton-antilepton 
pair, it would lead to decay of r-lepton into semileptonic channel, t" ^ l~+ mesons. The 
branching ratio (BR) for this inclusive process is expected at the same level as the BR for 
r ^ 3e decays, or somewhat smaller, because a decay of BH into a pseudoscalar state (e.g. 
TT-meson) may be suppressed and the dominant decay channel should be a scalar one. 

There can also be interesting modes of r decays with non-conservation of baryon and 
lepton numbers, as e.g. r~ — > e~e^p, e~e~p and analogous ones with neutrons and neu- 
trinos. The diagrams describing such decays are essentially the same as Fig. [T]i, with 
replacement of the initial fj, with r and with a BH which emits three quarks and one lepton 
instead of two leptons. The effective Lagrangian describing e.g. such decays has the form: 



where different V's are spinor wave functions/operator of the corresponding fermions and 
Qn is the coupling constant of BH to n fermions. According to the arguments presented in 
the preceding sections, creation of BHs in multiparticle collisions should be suppressed due 
to the necessity for several particles to meet in the same small volume and thus, e.g. (74 
should contain the BH radius to the fourth power, while the rest of the necessary dimension 
is supplied by the BH mass: 




BH 



(21) 



54 = R%Mbh, and §2 = RsMbh- 



(22) 
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The amplitude of the decay r pll is determined by the matrix element of this operator 
between an initial r and final 3Z + p states. The matrix element of making a proton out of 
three quarks is: 

TTl 

(piV'glvac) ~ J^Up, (23) 

where m.q ~ 300 MeV is the constituent quark mass, that is the characteristic energy scale 
of the process, and Up is the spinor wave function of the proton. 

We do not distinguish in this section between the masses of light constituent quarks, 
ruq ~ 300 MeV and the characteristic scale of strong interaction, Aqcd ~ 100 MeV. 
However, in future sections of the paper, where we consider proton decay and neutron- 
antineutron oscillations, the difference may be important. 

It is straightforward now to calculate the decay rate 

r,..o^.o-3(ii)^(.f)^(5)'(^)-^™. (-) 

where the first factor, lO^"', comes from three-body phase space. The life-time with respect 
to this decay would be extremely long: for M^, ~ 1 TcV and n = 2, we obtain r ~ 10^^ 
years. We have omitted here the loop counting factors but the effect is tiny even without 
them. 



5.5 i^T- mesons decays 

Similar considerations can be applied to rare decays of i^-mesons. In particular, if two 
quarks constituting i(''^-meson might form BH, this BH could decay into any neutral combi- 
nation of two leptons, e^e^ ji^ and /x^e^. The amplitude of this decay can be estimated 
as: 

A(/f^2/)«,.„(^)'^-^*, (25) 

where gxqq is the Yukawa coupling constant of JC-meson to ^g-pair and the third factor came 
from the integration of the loop of qq pair which combines into if-meson. The corresponding 
decay width is given by 

V{K' It) = '-k^ (HhV (26) 
^ ' 4(27r)5 \mJ \mJ ^ ' 

We have assumed that Mbh = and E = niK since no energy is taken away to other 
particles. With g^qq ~ 1, = 300 MeV, M* = 1 TeV and n = 2, we obtain for the 
life-time of this decay 

r(K° ^ II) ^ 0.16 s. (27) 

For K2 it gives the branching ratio BR(K^ II) ~ 10^^, which is 4-5 orders of magnitude 
above the experimental bounds: for example, BR{K'^ e+e") < 10~^^ and BR{K'^ 
e'^/j.^) < 5 ■ 10"-*^^. To avoid the contradiction we should either take M* > 3 TeV or to 
see if there is a possible suppression mechanism for this decay in this scenario. In fact, 
there is one. It is natural to expect that BH should have the quantum numbers of the 
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vacuum, i.e. it should be a scalar object. Hence the /C-meson, which is a pseudoscalar, 
cannot transform to BH directly, but should emit some other particle in such a way that 
the remaining combination of the quark- antiquark system would be scalar. The simplest 
way is to emit a vr'^-meson, while the remainder would make a BH which would decay into 
II. This mechanism means in particular that two body decays K ^ II are suppressed and 
not dangerous. 

The amplitude of the decay K nil is described by the diagram in Fig. [2] and is equal 

to 

<79 QKttS Tni - 

where gx-wSi the coupling constant of K and vr to the scalar state of quark- antiquark pair, 
has dimension of inverse mass. The factor m^/(27r)^ comes from the quark loop. The mass 
of BH is M'^^ = {pi^ +Pi2f' = rn\ + m1. — ImKET^. The life-time with respect to this decay, 
after integration over 3-body phase space with Mbh depending upon the pion energy, 
is equal to: 

t{K ^ nil) = 0.85 • 10^ s {gK.sm^r^ [^^^ J 



'300MeVy 6.4-10-3 
mq J fn ' 

where fn is related to integration over phase space: 



(29) 



/■(l+/^')/2 ^ .,_2_ 

fn= dx x/x2 - ^2 (l + ^2 _ 22.) +„+i (39) 

Here ^ = m-njmK- The factor 6.4 • 10-3//„ is equal to 1 for n = 2, to 0.82 for n = 3, and 
to 0.58 for n = 7. 

The experimental bounds on the branching ratios of — > ^^21 decays are about 
(3 - 5) • 10-1°, while its life-time is 5.2 • 10-*^ sec. Thus for n = 2, = 1 TeV, and 
Qk-kS'^-k ~ 1 we are on the verge of discovery of such decays. For larger n a larger M^^ or 
smaller gK-KS are needed. 

Similar estimates can be presented for decays of charged i('-mesons, e.g. to 

^ 7r+;u+e-, 7r+2i/, etc. (31) 

Experimental bounds on the branching ratios of these decays are about 3 • IQ-^^ [8]. The 
absolute probability of decay (j3ip is approximately the same as that of K'^ — > n^ll. Since 
the total life-time of is 4 times shorter than that of /C^, the predicted branching ratios 
are 4 times smaller and, according to the discussion above, they may be close to the existing 
bounds. 

This model has a few interesting features/signatures. The dominant anomalous decay 
mode is three body. The charge of the emitted pion is the same as the charge of the initial 
K. The probabilities of the decays with charged and neutral leptons in the final states 
are approximately the same. The rather large magnitude of the branching ratios of these 
anomalous decays of ii'-mesons make them very interesting/promising candidates in the 
search for non-conservation of global lepton quantum numbers. 
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Figure 1: a) Muon decay /x — > 3e with BH intermediate state, b) Diagram for the process 




Figure 2: Kaon decay K ti . i^'-meson probably cannot transform directly to a 

BH, because it is a pseudoscalar particle; instead, the emission of a vr'^ leaves a scalar qq 
system. 

Moreover, we note another appealing feature of all the decays with a BH intermediate 
state: if undiscovered weakly interactive light particles exist, such as for example possible 
sterile neutrinos or axions, they should be emitted by the BH with the same probability 
of the other light particles, if compatible with the BH quantum numbers. In fact, tiny 
BHs may provide a unique opportunity to discover such a kind of particles because, if the 
so-called "equivalence principle" holds, gravity couples to any form of energy with the same 
strength and does not distinguish one type of particle from another. Obviously, the emitted 
weakly interactive particles cannot be seen by the detector and the event appears as a 
process with violation of energy, momentum and, possibly, angular momentum. 



6 Baryon number violating processes mediated by black holes 



6.1 Proton decays 

Proton decay is described by a more complicated diagram, as shown in Fig. [3j To create an 
electrically neutral, colorless, nonrotating BH, one of the three quarks in the proton must 
emit a lepton pair through 7, W , or Z exchange; one of these leptons, together with the 
three quarks, may then form such a BH devoid of any quantum numbers. 

Since the BH is now formed in a 4-body collision, the probability of BH creation would 
be suppressed by an additional small ratio squared of BH volume to the proton volume. 
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The ratio of volumes appears in the second power because we have two more particles 
in the initial state in comparison with BH creation in /i — > 3e decay. In the latter case 
the BH could be produced in a two-body collision. This leads to the suppression factor 
(Rs-^-qcd)^ , where Aqcd ~ 100 MeV is the inverse proton size (in what follows we skip 
the sub-QCD). In fact, the situation is more complicated and a one parameter description 
may be impossible because the proton size is roughly the inverse pion mass, while the 
characteristic quark energies in the proton, which determine the coupling to BH are of the 
order of the constituent quark mass, i.e. ruq ~ 300 MeV. 

We will come to the same conclusion studying the transformation 3qe~ — > qq through 
a virtual BH, as is considered in Sec. [5j The amplitude of the reaction is proportional to 
92 gi/M^jj = Rg, where gn is the coupling constant of BH to n fermions; g2 is dimensionless, 
while 5(4 has dimension of inverse mass cubed, see Eq. ()22p . (Notice that in the case of 
/i — > 3e decay the amplitude is proportional to g2 and that's why the decay probability 
should be much larger.) 

The amplitude of the proton decay corresponding to Fig. [3] is equal to: 



Aip ^ l+qq) = - In ^ ^-f-^ i'li^pi'gi^g • (32) 
It leads to the decay rate: 



r(P-^^.-.) = ^(ln^) (^) [jfj I d.Al-2.r-^. (33) 

Correspondingly the lifetime of the proton with respect to the inclusive decay p qql^ is: 
.,.:0-.e«s(^)'°^^(2^)^-* (155M^)',„-.(„./XeV) ,34) 

where fp{n) is 1, 1.3 and 2.2 for n = 2,3 and 7 respectively. 

The best experimental lower bounds at the level Tp > 10^^ years f8] are established for 
the modes p — > e~^7r^ and p — > uK~^. For all other 2-body and some three-body modes the 
bounds are at the level of 10^^ years. The disagreement of our result with experiment can 
be easily avoided if we take a slightly larger M^,, still even smaller than 3 TeV. On the other 
hand, we should keep in mind that our estimates are by no means rigorous; they are only 
true up to an order of magnitude (and possibly some omitted factors would allow M^, to be 
as low as 1 TeV). For example, if we take into account that the proton inverse size, A ~ 100 
MeV is three times smaller than the quark mass rUg ~ 300 MeV, the life-time might become 
larger by the factor (m^/A)^ ~ 10'^. 

There is one more argument indicating that the decay life-time (j34p may be underesti- 
mated. By assumption, the intermediate virtual BH is a scalar and since we believe that 
gravity does not break parity, such BH cannot go into a pseudoscalar particle. In other 
words, the pair qq cannot form a pseudoscalar meson or vector meson. Thus we come to 
the important conclusion that proton must have predominantly 3-body decay modes: lepton 
plus two mesons with the two mesons in scalar state. This immediately shifts our estimates 
for the probability of decay into 2 mesons) down by an order of magnitude because 

of smaller phase space of scalar state. On the other hand the decay into three leptons is 
not influenced by this argument and its life time should be given by Eq. (I34p . Bearing 
in mind that the experimental lower bounds on the proton life-time with respect to three 
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BH 



Figure 3: Gravitationally induced proton decay. Since a 4-body collision is required in 
order to form a BH devoid of any quantum number, the process is strongly suppressed and 
experimental constraints can be compatible even with a gravity scale in the TeV range. 

body lepton channel, l^l^l^ (/ = e, /j,) are at the level (8 — 5) • 10^^ years we see that proton 
decays are on the verge of experimental discovery if is slightly larger than or about 2 



There are quite peculiar signatures specific to the model of proton decay considered here. 
First, as we have already mentioned, the decays should be mostly 3-body ones. Second, 
the final state particles must always contain a positron, e"^, or a positive muon, fi'^. The 
branching ratio into three lepton channel is predicted to be larger than that into e"*" and 
two mesons, because it is natural to expect that the probability of BH decay into a neutral 
combination of two leptons (or antileptons, or lepton and antilepton) is more or less the 
same as the probability of the decay into two quarks, while the probability of the subsequent 
quark transformation into two mesons is smaller than one, because other channels are open. 
The energy spectrum of emitted leptons (with the same charge as proton) is cut-off at higher 
energies due to the factor (1 - 2E/mp)^+^°/("+^), see Eq. (f33l) . 



6.2 Neutron-antineutron oscillations 

Another process where non-conservation of baryons is actively studied by experiments is 
neutron-antineutron transformation. While in many cases, as e.g. SU{5) GUT or elec- 
troweak theory, neutron-antineutron oscillations are impossible or completely negligible, 
because they demand change of baryonic number by 2, AB = 2, the gravitational breaking 
of global symmetries does not respect any selection rule and the oscillation time may be rea- 
sonably small. In the framework of the approach presented here, the neutron-antineutron 
oscillations are described by the diagram of Figure [H An estimate of this diagram is very 
uncertain and the result should be taken with great caution. There are two loops containing 
weak W 01 Z bosons. Both these loops are logarithmic and if the ultraviolet cutoff is given 
by the effective Planck scale, M^,, their contribution is not suppressed as an inverse power of 
the weak boson mass. Logarithmically divergent part of such loop diagram gives the factor: 



where Ajjv is the ultraviolet cutoff, which is probably reasonable to take equal to the 
effective Planck mass, Ajyy = M*. The amplitude of neutron-antineutron transformation 
contains this factor squared. 



TeV. 




(35) 
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The other part of the diagram, containing the lepton loop, is hnearly divergent after we 
perform integration in the loops containing weak bosons. The integral should be cut-off at 
the same energy scale as the neighboring triangle diagrams with weak bosons, i.e. at M^,. 
In other words, the linear divergent part as usually vanishes and the integral is proportional 
to the external momentum which, in this case, is of the order of M^. (See discussion in Sec. 
[5l two paragraphs below Eq. ([20])). 

As a result of this rather frivolous estimate we obtain for the amplitude of the transition 
of three quarks into three antiquarks: 



La 



B=2 



^In 



Ml 



(36) 



where the first factor comes from two triangle parts of the diagram, Eq. ()35p . the second 
factor is the ultraviolet cut-off of the linearly divergent loop with three lepton lines and 
one virtual BH and the next factor is the coupling constant of BH to 4 particles. The last 
term is the product of 6 quark operators and C is the charge conjugation matrix. 

Now it is straightforward to obtain the time of neutron-antineutron oscillations taking 
the matrix element (?^|i^A_B=2|^)• Since the effective energy cutoff in 3(7-transition into 
neutron is the QCD scale. A, we obtain: 



vr \mz 



(Ml 

V A 



7+; 



A" 



(37) 



(note that we have omitted here the huge loop counting factor (27r)^, because the result 
is weak anyhow). With n = 2 and M* ~ 1 TeV and A = 100 MeV it corresponds to an 
oscillation time of about 3 • 10^^ s which is twelve-thirteen orders of magnitude below the 
existing experimental limit: direct searches for n — > n processes using reactor neutrons put 
the upper limit on the mean time of transition in vacuum [8]: 



Tnn > 8.6 • 10^ s, 
while the limit found from nuclei stability is slightly stronger: 

Tnn > 1.3 • 10^ S. 



(38) 



(39) 



If the theoretical prediction of Eq. (j37p were true, the chances to observe (n — n)-oscillations 
in the reasonable future are negligible. However, one can obtain much more optimistic 
predictions if there exist supersymmetric partners of the usual particles, as considered in 
the following subsection. 



6.3 Supersymmetric extension 

Some of our estimates, such as for (n — n)-oscillations, would be different if there exist 
supersymmetric partners of the standard model particles. Since supersymmetry remains 
a hypothesis, not yet proven by experiment, and nothing is known about the masses of 
superpartners, except for lower limits on their values, more ambiguity is introduced into 
the calculations. Consequently we here consider anomalous processes with inclusion of 
superpartners as a separate subsection. 
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BH 

Figure 4: (n — n)-oscillation mediated by a virtual BH. If we consider only Standard Model 
particles, the effect is negligible. 

The spins of the SUSY particles may differ by 1/2 from their standard model partners, 
with all other quantum numbers being the same. In particular, there could be scalar quarks 
(s-quarks) or spin-1/2 partners of vector bosons mediating interactions. Existence of new 
types of elementary particles would modify both neutron-antineutron transformation and 
proton decay. 

In this case, one of the quarks in the neutron can emit a neutralino, x^j s-^id become 
a squark, q. This g, together with remaining quarks, can form a neutral and spinless BH. 
This BH in turn may decay into two antiquarks, 2g, and anti-squark, q. The latter captures 
X*^ and becomes the usual antiquark, q. This completes the transformation of three quarks 
into three antiquarks (see Fig. [5]). The analogous process with emission of a gluino does not 
help, because after this emission the remaining two quarks and one squark become colored. 

To find the amplitude corresponding to this diagram we need to calculate the contri- 
bution from the loop containing two s-quark propagators, one neutralino propagator, and, 
most problematic, a propagator of the virtual BH, about which we do not have much infor- 
mation. Possibly the result would be less ambiguous if we were to use the old non-covariant 
perturbation theory with particles on mass shell with positive energies. The last condition 
is important for definition of the vertex of interaction of BH with particles entering into 
it or emitted by it. Because of that, the total energy of incoming or outgoing particles is 
always of the order of the energy of the initial state and the mass of BH is of the same order. 
The necessity to use non-covariant perturbation theory in our description of interaction of 
particles with virtual BHs leads to breaking of Lorentz invariance. Another argument in 
favor of non-covariant perturbation theory is the impossibility to make Wick rotation with 
virtual BH - at least we do not know how to do that. 

The effective Lagrangian corresponding to the diagram in Fig. [5] is the following: 

where = R'^Mbh is the coupling constant of BH to 3 particles, two of which have spin 
1/2 and one is scalar, a ~ 0.01 is the electroweak coupling constant at characteristic SUSY 
scale; is the mass of neutralino, msusY is the mass of other superpartners, and we 
assume that they are all of the same order of magnitude, ~ n^susY- A subtle point is 
the value of the BH mass. According to the arguments presented above we take it to be of 
the same order of magnitude as the energies of the external particles; again, in the absence 
of any fundamental theory describing behavior of virtual BH, we can consider this at best 
a plausible assumption. 
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BH 

Figure 5: (n — n)-oscillation with supersymmetric particles. In this case the observation of 
the phenomenon may be accessible to future experiments. 

Taking the matrix element of operator (j40p between n and n states we find: 

{n\L^,^,\n) - - (2,^8 ^^^^^ [j^ ) • (^D 

Here A is the effective energy which enters in calculating of the matrix element 3 quark 
operators over the neutron state, it is usually taken about 100 MeV; and (27r)^ is the loop 
counting factor - there are 4 loops with virtual quarks which go either into n or n and each 
loop provides with l/(27r)^. 

Taking all the factors together we find for the time of (n — n)-oscillations: 

(5 4 4(n+2) 

r..«3.10 seclO".. (^^J (sOOCiv) J (tW J " ^^^^ 

This result looks quite promising. If is not too much larger than 1 TeV and the 
SUSY partners are not far from 300 GeV, the chances to observe neutron-antineutron 
transformations are very good. According to the model presented here their observation 
would indicate not low scale gravity but also low energy SUSY; but this is probably too far 
fetched. 

We also note that the contribution of SUSY partners to proton decay is negligible. 



7 Heavy quark decays 

We discuss here decays of heavier quarks or mesons containing such quarks. Though the 
experimental accuracy in their decays are much lower than for proton or muon, the effects 
may be amplified because of larger masses and it is quite probable that the manifestation 
of low scale gravity will be first observed in such decays. 

Let us start from the heaviest, the t-quark. Since it has a very large mass, mt = 175—180 
GeV, close to the assumed gravity scale, one can expect that B-nonconserving decays of the 
t-quark may be noticeably enhanced. A B-nonconserving decay t — > 4q + / is described by 
a diagram of the type presented in Fig. [6l The same considerations as those presented in 
Sec. 5 lead to the decay rate: 
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where £5 ^ 10 is the 5-body phase space normalized to the t-quark mass. 
The total decay rate of the i-quark is known to be: 

Ttot ^anit^r- W-^mt. (44) 

Using Eqs. (|43|) and (j44|) we obtain the branching ratio (always assuming M^, ~ 1 TeV): 

BRaBj^o ~ 10~^° - 10"^^ (45) 

for n between 2 and 7. Probably the decay rate (jl3]) is underestimated because of too large 
phase space suppression factor and the branching ratio is somewhat larger. In particular 
it would be larger for B-nonconserving decays of mesons or baryons containing t-quark 
because phase space suppression in this case would be much milder. On the other hand, 
the loop counting factor, omitted above, would play its destructive role. 

At the present time there are no experimental restrictions on nonconservation of the 
baryonic number in decays of particles containing t-quark. However, a noticeable baryonic 
charge non-conservation in t-quark decays and some other anomalously enhanced decays of 
mesons containing t-quark may be accessible to future experiments. 

More realistic possibility may be in the leptonic sector: for example, the BH in Fig. [6] 
can emit an e^/i^-pair, leading to the violating the family lepton number decay t ue/i. 
Its amplitude is 

Ait ^ ueti) = In RliPuAi^^^e. (46) 

The predicted decay rate is 



r(*-».) = ^f^>n^(^j (^J -I (47) 

and, for n = 2, we expect a BR at the level of 10~^. 

On the other hand, t-quark may be not the best candidate for search of anomalous 
decays induced by virtual BH because rrit > mw and the total decay width of t-containing 
particles is quite large due to the open channel into W and 6-quark. 

Probably a better place to search for low scale gravity effects could be decays of 6-quark 
or, to be more precise, decays of mesons containing 6-quark, in particular i?-mesons. On 
one hand, 6-quark is lighter than t-quark, its mass is about 4.5 GeV, and this makes the 
effects weaker. On the other hand, the total decay widths of particles containing 6-quark 
are smaller than those with t-quark (because h is not heavy enough to decay into W boson 
and lighter quarks) and this would enhance branching ratio of anomalous decays induced 
by gravity. 

As an example, let us consider a decay of iJ'^-meson consisting of a heavy 6-quark and 
light d-quark. As is known from QCD such a system of light and heavy quark has the size 
of the order of A~^ and the characteristic energy of the light quark about A. Its decay into 
two light leptons can be considered in the same way as muon or IT- meson decays in Sec. O 
One should expect that an uncharged and nonrotating virtual BH can be formed directly 
in collision of b and d and it is not necessary to invoke any other virtual particles. If this 
is true (but remember that BH should probably be a scalar and not a pseudoscalar), then 
the amplitude of 5-meson decay into two quarks or two leptons is given by the expression 

..R „N aloBqq 4 [ d^p (ppi + mbmq) - 

-"t^-")- 147(2^/ C**) 
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where gsqq is the couphng constant of i?-meson to two quarks and pi = p — Pb- 

The cut-off of this integral is determined by the strong interaction scale A and could be 
described by a formfactor: 

Fi = F{pl/ml,pl/ml,pl/ml) (49) 

which vanishes if the participating particles are too far from the mass shell. However, it is 
rather complicated to impose this condition in the standard form of the Feynman integral 
We will use to this end the on-mass-shell representation of the Green's functions: 

G{^,t) ^ J j^^exp{-iEt + ipx), (50) 



where E = \/ p'^ + m?. The one loop diagram is now described by the expression 

PBVb + rribmg) . (51) 



Eq Ef, {Eb — Eg — Eij) 



The last term in brackets comes from the fermionic trace and contains the product of four- 
momenta of the virtual fermions. The new formfactor F2 is a function of three-momentum 

because all the particles are on-mass-shell but energy is not conserved in each vertex. 
This form-factor is determined by the interaction potential and is cut-off at |p| ~ A. 

This one-loop integral (15ip can be easily estimated giving the result ~ A^. An important 
fact is that it does not contain the mass of heavy quark, in the denominator. 

The decay width of i?-meson into the channel B ^ II 01 B ^ qq can now be estimated 

as 

For gsqq = 1, n = 2, A = 100 MeV and = 1 TeV we find the hfe-time ~ 3 • 10"^ 
s. The total life-time of B^ is 1.5 • 10"^^ s. Thus the branching ratio of anomalous decays 
with the chosen values of the parameters should be about 5 • 10"^*^. This result is below 
the existing experimental bounds. The branching ratios of B^ into e^e~ , ^jl^ ^~ , efi are all 
bounded by approximately 10"''. We repeat, however, that estimated branching ratio may 
be true if pseudoscalar BHs are allowed. 

If the virtual BH must be scalar, and we consider the decay B 21, then before 
collapsing into BH the system bd should emit a light pseudoscalar (PS) meson, P, which is 
later to be absorbed by the final state leptons, but the probability of that is very low (it is 
a weak interaction process). Thus the 3-body decays should dominate in the same way as 
found for K-meson decays in Sec. [5l 

The amplitude of the decay B P + 21, where P is a light PS-meson, is determined 
by the similar loop integral as above with the only change that m?^ — 2Epmh is substituted 
instead of m^. This leads to further deviation from the mass-shell pole and the integral is 
suppressed by approximately an order of magnitude in comparison with the previous case, 
i.e. ~ O.lA^. Correspondingly the decay width is equal to 



r(...o. '°----(— '- (A)- (-)-^ 



(53) 
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Figure 6: Top decay with BH intermediate state and violation of baryon and lepton numbers. 

where qbps is the couphng constant of transition of i?-meson into PS-meson P and a scalar 
state of two quarks. It has dimension of inverse mass and it may be natural to assume that 
9bps ~ 1/^ (this is a quite strong coupling). 

For gsps ^ l/K n = 2, K = 100 MeV and = 1 TeV the hfe-time with respect to 
this decay is 0.2 s. It leads to the branching ratio BR{B — > Pll) ~ 10~^^. This relatively 
large branching is related to a huge coupling qbps ~ 1/^- 

An interesting process could be a decay of into p + e~ . The branching ratio of this 
decay should be suppressed with respect to the decay of B^ into two leptons ([52]) at least 
by the factor (i?5A)^ < 10~^^, because now the BH emits four particles and therefore the 
amplitude of the process is proportional to 9294,/ M'bh- makes this decay impossible to 
observe in foreseeable future. 

There can be some other i?-meson decays as well, with violation of baryon or lepton 
numbers or not, where virtual BH could give noticeable contributions and an observation 
of possible anomalously large branching ratios might indicate on the effects of virtual BHs. 



8 Conclusions 

Observable effects of gravity in particle physics are an interesting and fascinating possibility 
allowed by models with large extra dimensions, where the fundamental gravity scale can 
be set in the TeV range. In such a framework, we have considered rare decays with BH 
intermediate states. Since BH decays violate global symmetries, we expect lepton and 
baryon numbers non-conserving processes. We elevate the classical requirement of zero 
charge and zero angular momentum for sub-Planck- mass black holes to a general conjecture 
even in quantum gravity. In that case, the predicted B or L non-conserving decay rates 
or the frequency of neutron-antineutron oscillations are interestingly close to the existing 
experimental bounds but not in conflict with them for M^, > 1 TeV . It is striking that 
neutron-antineutron oscillations and some anomalous decays of K and i?-mesons are very 
close to the existing experimental bounds and, thus, quite promising for observation of TeV 
gravity effects. 

We have taken M^, as low as possible, namely about 1 TeV. With higher M* all exper- 
imental bounds can be easily satisfied, but new physics becomes harder to test. Unfortu- 
nately we do not know how large can be M^,; hence it is impossible to reject our model by 
experiment. Even with slightly larger than 1 TeV all effects discussed here would be 
strongly suppressed. 
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The model presented here is very speculative and includes plenty of wishful thinking. On 
the other hand, it gives very interesting testable predictions with specific signatures absent 
in other models. We repeat that with the "natural" value of = 1 TeV the predictions 
for the magnitude of new effects are quite close to the existing experimental accuracy. 

The way in which virtual BHs are treated here implies some striking features, e.g., 
possible breaking of Lorentz invariance. It may mean, in particular, that the magnitude of 
effects in different coordinate frames may be significantly different - a reincarnation of the 
old concept of the ether. 
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Appendix: Short review on current constraints on 

To begin with, we stress that, here and throughout the paper, as fundamental gravity scale 
we take the quantity M*, which is related to the (4+n)-dimensional gravitational constant 
of the (4+n)-dimensional Einstein-Hilbert action by the simple relation 

M:+2 = 1 . (54) 

In the trivial case n = 0, we have = Mp^ = 1/Gn- In the literature there are at least 
other two popular conventions. Some details on the different possibilities can be found, for 
example, in Ref. [29j. 

In addition to this, we now assume that the (4+n)-spacetime is given by x T", 
where M4, is the standard 4-dimensional spacetime we know and T" is an n-dimensional 
torus of radius R. In this special case, the volume of the extra dimensions is finite and 
equal to (27ri?)"' and the relation between the standard 4D Planck mass Mpi = 1.22 • 10^^ 
GeV and the fundamental gravity scale M^, is 

Mil = (2^^)" ^r"^^- (55) 

With these two statements in mind, we now review present experimental lower bounds on 
the magnitude of M*. Prom now on, we follow Ref. [8]. 

Since a TeV gravity scale allows the emission of gravitons at colliders [39j . constraints 
on M=K can be obtained looking for missing energy in processes such as e^e~ 7G (the 
probability that the graviton G interacts with the detector is suppressed by the standard 
Planck mass Mp;, and therefore negligible). The non-observation of such events at LEP 
leads to the 95% CL bounds 00] 

> 1.43, 0.76, 0.47, 0.33, 0.25 TeV (56) 

for n = 2 — 6. 

Much more stringent constraints for n < 4 can be obtained by astrophysical consider- 
ations. We note however that these bounds require the existence of gravitons lighter than 
about 100 MeV and can be evaded if gravitons acquire small extra contributions to their 
masses, as suggested in [HI 321 US] (in the scenario of Ref. jl3] , even the n = 1 case can- 
not be safely excluded). In fact, the astrophysical environments used to constraint M* are 
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characterized by a typical energy per particle of 10 - 100 MeV and an effective graviton 
mass of 100 MeV or more prevents a copious production of them. These scenarios leave 
instead unchanged collider constraints, where the typical energy is roughly 100 GeV and a 
100 MeV mass is unimportant. 

Since gravitons are weakly interactive particles, a possible their production in a super- 
nova event can compete with neutrino cooling. Neutrinos detection from SN1987A requires 

m 

> 12.5, 1.0 TeV (57) 

for n = 2 and 3 respectively. 

On the other hand, if we consider all gravitons produced by all the supernovae in the 
history of the universe, we can expect a diffuse gamma ray background due to the graviton 
decay into photons. The non-observation of such a diffuse background by EGRET satellite 
puts the bound [l5] 

> 34, 2.6 TeV (58) 

for the n = 2 and the n = 3 case. 

Finally, noting that gravitons produced in supernovae events are not high relativistic 
particles, we can expect that many of them remain gravitationally bound to the neutron 
star remnant and that their subsequent decay reheats the surface of the star. The measured 
luminosity of some pulsars leads to the very stringent bounds |46] 

> 670, 20 TeV (59) 

always for n = 2 and 3. 
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